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1. Introduction 


The radiation dose received from high energy galactic cosmic rays (GCR) is a limiting 
factor in the design of long duration space flights and the building of lunar and martians 
habitats. It is of vital importance to have an accurate understanding of the interactions of 
GCR in order to assess the radiation environment that astronauts will be exposed to. 

Most previous studies have concentrated on strong interaction processes in GCR. However 
there are also very large effects due to electromagnetic (EM) interactions. EM studies have 
previously concentrated on single photon exchange leading to nucleon removal. However two- 
photon processes also occur which lead to the production of lepton pairs with cross sections 
of the order of kilobarvs. Also at high energy the stopping powers from these processes can 
exceed that due to atomic collisions. Thus even though very high energy GCR are not as 
abundant as lower energy GCR they still must be considered due to the fact that the cross 
sections and stopping powers are so much larger than normal. 

In this report we describe our first efforts at understanding these EM production processes 
due to two-photon collisions. More specifically, we shall consider particle production 
processes in relativistic heavy ion collisions through two-photon exchange. Examples of 
this broad category of processes include: 


z x z 2 — ZxZ 2 i + r 

(1.1a) 

Z\Z 2 — ► Z\Z 2 s + s~ 

(1.1b) 

ZiZ 2 — ZiZ 2 V + V~ 

(1.1c) 

Z Y Z 2 -> ZiZ 2 tf° 

(l.ld) 


in which l + l~ denote charged leptons, s + s~ denote charged scalars, V + V~ denote charged 
vector particles, and H ® is a neutral Higgs scalar. 

We shall limit our consideration to cases in which the colliding nuclei are identical, so 
that Z\ = Z 2 = Z. An important Feynman diagram that contributes to (1.1a), (1.1b), and 
(1.1c) is shown in the following figure (fig. 1). 



For process (l.ld), an important diagram is shown in figure 2, in which the triangular 


loop receives contributions from quarks, leptons and W gauge bosons. 

These processes are important for the following reasons (ref. 1). 

(1) These kinds of processes become increasingly important as the energy of the colliding 
nuclei increases, since their cross sections increase with energy. Thus their contributions 
to the stopping power of high energy ions also become more important at high energies. 

(2) These processes can be channels for production of charged particles, e.g., l + l~ , W + W~ , 
and neutral particles such as Higgs bosons, and various mesons. 

(3) For high Z nuclei these processes can be used for studying non-perturbative effects in the 
electromagnetic interaction. 

(4) They must be taken into account in the study of strong interaction effects in heavy ion 
collisions since they can lead to important background events, and must be taken into 
account also in the design of experimental set up, since they can lead to significant beam 
loss. 

Section 2 of this report gives a brief survey of a few major approaches used in the 
calculations for these processes. Section 3 examines some results of our calculations. We 
then point out briefly some open questions and make a few concluding remarks in Section 4. 

The purpose of this report is threefold. (1) It gives a simple, elementary introduction to 
this field. (2) It provides sample calculations for illustrating the approach we use. (3) The 
background and techniques developed here can be used as a general base for launching further 
and more specialized studies into this field. 

While it is not our main goal here to obtain new and original results, some of our results 
are possibly new, and are as yet not available in the literature. 


2. A Brief Survey of Different Approaches 

In this section, we briefly list a few major approaches used in calculating cross sections 
for the kind of processes we are interested in. The first approach has been discussed in 
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references 2 and 3. In this approach, each colliding nucleus is replaced by an equivalent 
spectrum of photons. Each nucleus is considered to move in a straight line, unperturbed by 
the interaction. At a distance b from the line of motion of a nucleus, a spectrum of photons 
is generated, whose frequency distribution has the form: 


N(u, b ) 



(2.1a) 


where 



(2.1b) 


Kq, K\ are modified Bessel functions, see reference 4, Sections 3.7 and 15.4. 


The cross section for this process can be written as an integral of a photon distribution 


function multiplied by a photon-photon cross section. 


a = 



du;2 

F(u)l, U>2 )°77 (u>i , ^ 2 ) , 

UJ 2 


(2.3a) 


where 


F(u>i,U 2) =2ir b\db\ 

~ ^lmin 



d<f>N{u)\, b\) 


x N{u 2 MW -R1-R2) 


(2.2b) 


and 



(2.2c) 


where lo\ and u >2 are the frequencies of the photons emitted by the nuclei, 61 and 62 are the 
distances of the nuclei from the point where the photons collide. Details can be found in 
Appendix A. Various differential cross sections can be derived from these equations. First 
we consider where W is the mass of the produced charged particle pair. We note that 
W 2 = Hence we can equate in (2.2a) 


a = 



(2.3a) 
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and 


da 
dW 2 


1 f du>\ ( W 2 \ ( W 2 \ 

~W^j Ui F ( U ’ 1, 4wJ <Tl ' 1 r 1, 4w 1 j- 


w 2 \ 


(2.3b) 


Next we define the probability for producing a particle pair P(b ) at impact parameter b by 

1 da 


p(b) = 


2irb db 


(2.4a) 


where 

^ = / — / ^F(v h v 2 )'7-,-,(“l,‘S2)HI> ~ b'), (2.4b) 

ab J uj\ J u 2 

in which it is understood that the 6 function is to be taken inside the triple integral which 
defines u^). The correctness of (2.4) can be checked by integrating both sides of (2.4b) 
over all values of the impact parameter b, which then yields (2.3a) for the total cross section. 
P(b) is the probability for the events in which two nuclei collide with each other at impact 
parameter b, producing a charged particle pair in the process. A quantity L, known as the 
two-photon luminosity function is defined by (see ref. 2, eqs. (1), (9), and (10)) 


So 



(2.5a) 


(2.5b) 


and 

da dL > _ _ _o . j . 

= dfV^^ (W 1 (25C) 

where we have used the fact that cr 77 ( u\,U 2 ) actually depends only on W 2 so that we can 
write 


a~n 


(u>i,W2) = o^iW 2 ). 


( 2 . 6 ) 


It is our view that equation (10) of reference 2 is in error, and have duly corrected the error 
in the above definition of the luminosity function L. For the stopping power calculation, we 
use the formula 


dE 

dx 



— ~{u\ +a;2)F(a;i,a;2)u'^(u;i,C4^), 

u>2 


(2.7) 
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where p is the number of nuclei per unit volume. 

The second type of approach has been applied to a related set of purely quantum 
electrodynamic (QED) processes: e + e” — > e + e~l + l~. This process can be calculated 
within the framework of QED. Cross sections can be obtained numerically by Monte-Carlo 
integration. Approximate formulas for total cross sections have also been obtained. See 
references 5 and 6. This kind of approach can be modified to apply to relativistic nucleus- 
nucleus collisions, provided one takes into account properly the effects of nuclear currents. 
See reference 7, Section II. 

In an approach closely related to this second type of approach, Bottcher and Strayer 
treated the colliding nuclei classically, by regarding them as classical charge distributions. 
The remaining amplitude for the production of charged particle pair is then obtained in the 
framework of QED. Thus for the case of the reaction Z\Z<i — *• Z\Z 2 l + l~ the total cross 
section can be written in the form (ref. 8, eq. (10), p. 38): 


a = 




u (p-> s -) 

t a . I- 


s_,s+ 


+ h 


i>-~ h ~ ™-i 


h 


i>-~ - "if 

)2 

v(p+,s+) 


h 


( 2 . 8 ) 


where v denotes the velocity of one of the nuclei in the center of momentum frame, p~ and 
p_l_ are the momenta of the produced leptons, s_ and s+ are their polarizations, u(p_, s_) 
and v(p+, s-|-) are the lepton spinors, k\ and k 2 are the momenta of the exchanged photons, 
and fi and /2 are the nuclear form factors. For any 4-vector A, the slash notation 4 is 
defined by 

4 = E VA < 2 - 9 > 

where 7^ are the Dirac 7-matrices. See for instance reference 9, Appendix 2, pages 355-361. 
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3. Results 


In this report, we adopt the approach discussed in references 2 and 3. As samples of our 
calculations, we present a number of results for the process 208 P 6 208 P 6 — + 208 P 6 208 P 6 l + l~ , 
and some others. Most of our calculations are done for colliding beam energies of 3400 Gev 
and 8000 Gev per nucleon. The impact parameter b varies over the range from 10 fm to 
1000 fm. The mass of l + l~ varies from a threshold equal to 2 ra/ up to about 1000 Gev. In 
Appendix B we list the photon-photon cross sections for the following processes: 


77 - l + l~ (3.1a) 

77 — ♦ s + s~ (3.1b) 

77 - V'+V'" ( 3 . 1 c) 

77 — * ► H° (3. Id) 


The derivations of some of these cross sections are also given there. By using ( 2 . 2 )- 
( 2 . 6 ), we can then obtain various luminosity functions, differential and total cross sections, 
probabilities, and stopping powers. 

Table 1 shows the total cross section for Z Z — ► ZZe + e~ . We compare our numerical 
results based on (2.1) and (2.2), with the results based on the FLacah formula (ref. 5, eq. (F.l), 
p. 276) 


a 


28 (ZiZ 2 a 2 ) 2 
277rm 2 


(, l 3 - Al 2 + BI + C) 


(3.2) 


where 


A « 6.36, B « 15.7, C « -13.8, 




Z{, pi, and m,-, i = 1, 2 are the charges, momenta, and masses of the colliding nuclei. 
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Table 1 


Colliding 

nuclei 

Z 

Incident 
energy /nucleon 
E (GeV) 

Total cross section (kilobarn) 

Our results 

Calculated 
from formula 

0-16 


0.2241 x 10" 1 

0.2216 x 10 _1 



0.3052 x 10 _1 

0.3020 x 10 _1 

Al-27 

3400.0 

■BHI 



8000.0 

■H 

B 

Fe-56 

3400.0 

0.2500 x 10 1 

0.2473 x 10 1 


8000.0 

0.3404 x 10 1 

0.3370 x 10 1 

Pb-208 

1 

0.2473 x 10 3 

0.2445 x 10 3 


| 

0.3366 x 10 3 

0.3333 x 10 3 


Table 2 shows the corresponding stopping power calculations. The energies of the incident 
particles are given for both the case of colliding beams and also the case of an incident beam 
colliding with a fixed target. 


Table 2 


Incident energy/nucleon 
E (Gev, colliding beams) 

Incident energy /nucleon 
E (Gev, fixed target) 

l v / dE\ 
P x 

(Gev fm 2 ) 

(~3f) 

(Mev/cm) 

0.9636 
0.1367 x 10 1 
0.2704 x 10 1 

3400.0 

8000.0 

1.039 

3.036 

14.64 

0.2462 x 10 8 
0.1363 x 10 9 

0.2129 x 10 -2 
0.4585 x 10 2 
0.6909 x 10 3 
0.2032 x 10 11 
0.1182 x 10 12 

7.02 xlO -4 
1.51 xlO 1 
2.28 xlO 2 
6.70 xlO 9 
3.90 xlO 10 

For Pb-208, p « 3.30 x 10 22 cm -3 
























Note that the stopping power for e + e" production increases with energy. So as the 
energy of the colliding nuclei increases, the contribution of this process to stopping power 
also becomes more important. In contrast, the contribution to stopping power from atomic 
collision and other processes first increases with energy, and then decreases. Hence, as energy 
increases, eventually these other contributions become less important. To put our results 
into perspective, we note that for Fe-56 at a kinetic energy of 1 Gev, its stopping power in 
water due to atomic collision is around 10 4 Mev/cm, see reference 10, figure 2.15, p. 74. 

In figure 3, we give plots of as a function of W in different ranges of W. The 

differential cross section can be obtained from by multiplying by a 77 cross 
section as in (2.5c). 

Figures 4a-d show plots of oyy(W 2 ) for the reactions 77 — ► l + l~, 77 — ► s + s~, 
77 _» v + V~, and 77 — ► H°. 

Figure 5 shows plots of P(b ) for the reaction 208 P 6 2 ^ 8 P& 208p^208p^ e + e - a t different 

energies. 

Figure 6 presents plots of the total cross section for the process 208 P 6 208 P 6 — ► 

208p 6 208p 6 H 0 

We have compared some of our results with the published results of Papageorgiu and 
Baur, and found good agreement. In the following, we give a sample of such comparisons. 
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Table 3 



W (Gev) 

w/^s 

tx /2 dL 

K dW? 

Our result 

Papageorgiu’s 

result 

Incident 

energy/ 

nucleon 

(colliding 

beams) 

E = 3400.0 
Gev 

100.0 

0.7070 x 10 -4 

0.3152 x 10 3 

0.33 x 10 3 

141.4 

0.1000 x 10" 4 

0.8630 x 10 2 

0.90 x 10 2 

212.2 

0.1500 x 10" 3 

0.1206 x 10 2 

0.13 x 10 2 

282.9 

0.2000 x 10“ 3 

0.1990 x 10 1 

0.21 x 10 1 

E = 8000.0 
Gev 

19.0 

0.5709 x 10" 4 

0.6129 x 10 3 

0.70 x 10 3 

280.0 

0.8413 x 10" 4 

0.1708 x 10 3 

0.20 x 10 3 

370.0 

0.1112 x 10 -3 

0.5444 x 10 2 

0.60 x 10 2 

460.0 

0.1382 x 10" 3 

0.1881 x 10 2 

0.20 x 10 2 

550.0 



0.70 x 10 1 

640.0 

0.1923 x 10“ 3 

0.2606 x 10 1 

0.28 x 10 1 
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Table 4 


Incident energy/ nucleon (colliding beams) E = 3755.6 Gev 

W(Gev) 

dw? (fm 2 Gev 2 ) 

Our result 

Baur’s result 

0.1200 xlO" 2 

0.6188 xlO 13 

0.62 xlO 13 


0.4923 xlO 13 

0.48 xlO 13 

0.1732 xlO" 2 

0.2777 xlO 13 

0.28 xlO 13 

0.2000 xlO" 2 

0.1714 xlO 13 

0.17 xlO 13 

0.2200 xlO" 2 

0.1222 xlO 13 

0.12 xlO 13 

0.2400 xlO" 2 

0.8894 xlO 12 

0.92 xlO 12 


0.6607 xlO 12 

0.69 xlO 12 

0.3000 xlO" 2 

0.3846 xlO 12 

0.40 xlO 12 


Papageorgiu and Baur’s results were taken from appropriate graphs in their papers (ref. 
2, fig. 3; and ref. 3, fig. 9). 

Cross sections are expected to scale roughly as Z\ Z\. For our case Z\ = Z 2 = Z. So in 
order to obtain the corresponding cross sections, luminosity function, or stopping power for 
different nuclei, one can simply multiply the results we have here by a factor ^—82)^ ' Thus 
if one wants the results for Al Fe collision, one can multiply the results presented in this 
section by the conversion factor The different nuclear sizes are expected to affect 

the results also. However for a rough order of magnitude estimate, such a simple scaling is 
expected to be reasonably accurate. 
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4. Open Questions and Conclusions 


For small values of b, and mj, such as m/ = m e , P{b) exceeds 1. This signifies the 
breakdown of perturbation theory. The question as to how to extract meaningful results 
from theory is under active investigation. See reference 11. In our simple approach, 
we have regarded the nuclei as point charges. By using form factors for the nuclei, the 
problem of violation of unitarity is expected to be somewhat ameliorated. However this 
problem still needs to be addressed, because for high Z nuclei, the coupling constant for 
the electromagnetic interaction is of the order Ze, even with nuclear form factor taken into 
account, which may therefore still lead to a breakdown of the perturbative approach to cross- 
section calculation. In a collaboration with Mirek Fatyga of Brookhaven National Laboratory 
(BNL), we shall investigate lepton pair production and neutral meson production (such as 
e°, 77 0 ) in high energy heavy ion collisions. In these processes, we shall look for possible 
deviation in the measured rates or cross sections from values calculated by perturbation 
theory. 

In many studies of the type of processes considered here, various approximations are 
used. We have mentioned the equivalent photon approximation, and the semi-classical 
approximation. Also, in the approach of references 2 and 3 , which we have adopted in 
this report, the effect due to phase coherence of the electromagnetic field generated by the 
nuclei has not been properly taken into account. One needs to investigate how valid these 
approximations are and what the regions of validity are for them. 

When one is primarily interested in the kind of electromagnetic processes discussed 
here, one needs to be able to estimate reliably the background due to strong interaction. 
Furthermore, there are other electromagnetic processes that also need to be studied, in 
addition to the ones we have looked at, even though the ones we have considered are among 
the most important. 

In summary, we have given a brief introduction to two-photon exchange processes in 
high energy heavy ion collisions. Our calculations are based on an approach discussed in 
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references 2 and 3. In view of the significance of this class of processes, and the many open 
questions that remain to be answered, we believe that further study in these areas will be 
valuable, not only for gaining a better understanding into these processes themselves, but 
also for studies and experiments in strong interaction physics. 

In the following Appendices, we discuss the derivation of some of the formulas we have 
used. We look at the equivalent photon approximation in Appendix A and show how this is 
applied to the two-photon exchange processes in relativistic nucleus-nucleus collisions. Then 
in Appendix B, derivations are given for some 77 cross sections. Appendix C provides a 
derivation of the fermion contribution to the process H ® — » 77 . In Appendix D, we look at 
the details of how certain integrals encountered in our calculations are evaluated. Finally 
Appendix E gives a simple derivation of the formula (2.7) used for calculating stopping 
power. 
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Appendix A. Equivalent Photon Approximation 


Consider a charge q moving along the x-axis. The effect of this charge on another charge 
located a distance b from the x-axis can be approximately calculated as follows. 

By first considering the electromagnetic (EM) field due to q in its own rest frame, and 
then making a Lorentz transformation to the laboratory frame, it is straight forward to show 


that the electromagnetic field due to q is given by 

Ei = -qv'ytib 2 + 7 2 u 2 f 2 ) -3 / 2 (A.la) 

E 2 = qby{b 2 + 7 2 u 2 t 2 ) -3 / 2 (A.lb) 

B 3 = -E 2 = q~ tn(b 2 + 7 Vt 2 )- 3 / 2 (A.lc) 

c c 

E 3 = Bi = B 2 = 0 (A.ld) 


t = 0 corresponds to the instant when q passes through the origin. When v as c, the 
components E 2 and £3 can be thought of as the components of a pulse of plane-polarized 
EM wave travelling along x. The energy flux of this EM field is given by the Poynting vector 


S = —E x B. 
47 r 


(A.2) 


So ignoring E\ for the moment, S points along x, and its magnitude is 


c *?2 


' s ' “£*• 


(A.3) 


in which we have made the approximation ^ ~ 1. Over a unit area, the flow of energy is 

/ oo r roo 

|S| dt = -j- I E 2 (t) dt. (A. 4a) 

-CO 47T J — 00 

Using Parseval’s theorem, we therefore have 

/ oo r roo 

\$\dt = / \E 2 {u)\ 2 dLJ, (A. 4b) 

-OO ^7T J — oo 

where E<i is the Fourier transform (FT) of i?2i defined by 


1 r °° 

s VK L m 


.-tut 


dt 


(A. 5) 
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Hence the quantity £2 (u>), defined by 


s 2 M = — iejMP, 


(A. 6 ) 


can be thought of as the energy per unit frequency per unit area of the EM field at frequency 
oj generated by the moving charge q. To obtain the photon number per unit frequency per 
unit area at frequency u>, we set 712(0;) = since each photon has energy hu. For 

the function 712(0;), the dependence on the distance b is implicit. To make the dependence 
on b explicit, we can write instead 

n 2 {u,b) = -^-S 2 {u). (A. 7 ) 


FVom (A. lb), we obtain 


Hence 



(A. 8 ) 


(A. 9 ) 


The remaining component E \ , of the EM field can be complemented by a magnetic field so 
that they can be considered to form a pulse of plane polarized EM wave. The same treatment 
can be applied to these components, so that the energy spectrum can be similarly obtained 
as before. The result is 


c , - 


SiM = -l-EiMf = 


-i-L-L 

47T 7 ^ i /2 


m 2 2 M\ 2 

\bJ 7 r / 





(A. 10 ) 


The effect of this pulse is roughly ^ that of the first pulse. So at high velocity, the second 
pulse can be neglected when compared with the first pulse. 

In conventional treatment, the two pulses are simply added together, so that the effect 
due to the original moving charge q is replaced by a spectrum of photons whose number 
density is simply the sum of the number densities from the two pulses discussed above. Thus 
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one set 


. v 1 r. i 

b) = + S\(uj) 


huAn v 2 


After identifying q = Ze, 


x 

e 2 

Tic 



= a, and noting 




= 5j(cj) for i = l,2, 


the photon energy spectrum 


(A. 11) 


(A.12) 


N(uj, b ) = + -^ 2 ( 0 ;) + 5i(— u) + S2(—v) 

= , 2 (u_\ 2 l 2 fwb\ 1 2 fub\] 

7r^ c \v ) \7v / [ * \7v / 7^ 0 \7V / 

- » ^ (- c ) 2 ( ^ ) 2 k ( -) + 4 W -) 

7T Z \v/ \ 7 V / L \ 1 V J \1V J 


(A- 13 ) 


Application of the Equivalent Photon Approximation to Two-Photon 
Exchange Processes 

When two nuclei Z\ and Z2 collide with each other, their EM interactions can be studied 
in terms of the EM interaction of the spectra of photons emitted by the nuclei. The situation 
can be pictured as in figure A. 2. 

The two photons 71 and 72 are considered as colliding head-on with each other. Taking 
a cross-sectional view perpendicular to the direction of motion of the nuclei, the situation 
can be pictured as shown in figure A. 3. 

Prom our previous discussion, the number of photons emitted by Z\ at P , whose 
frequencies are between u\ and u>\ + du\, is n(u)\,bi)duj\ 61 db\ d<f>\, where n{u)\, b\) 
is defined by (A. 11). Similarly, the number of photons incident at P emitted by Z2 is 
n(uj2,b2)duj2 f>2 <#>2 ^ 02 - Therefore the EM cross section for the collision of Z\ and Z2 
through two-photon exchange can be written as 


<7 = J n( ul M a-y-y(u!i,uJ 2 )bi db\ d®\ 62 df>2 d 4 > 2 

x 9 (b — R\ — R.2)dbj\ du2 


(A. 14 ) 
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in which 72 1 and i?2 stand for the nuclear radii of Z\ and Z 2 , and the 0- function takes into 
account that when b < R\ + R 2 , the two nuclei overlap, and the EM interaction is swamped 
by the strong interaction of the nuclei, and so one needs to restrict b to values > R\ + R 2 if 
one wants to look only at EM interactions. 

Since 6 = + bj - 2bib 2 cos 0) 1 / 2 , the integration / d(f> 1 d<f > 2 in (A. 14) can be simplified 

if one integrates over <p\ and converts the integration over <j > 2 into an integration over 4>: 


J d<t>\ d<f>2 


— ► 2n 


/ 


d<j> 


(A.15) 


So (A.14) can be rewritten as 


(A. 16) 


c — 2 tt J n(uJi,bi) n(u 2, 62) c r 77( w li a; 2)^(f > — R\ — R%) 
x 61 db\ b 2 df>2 d<p du>i du 2 - 

If one now substitutes for n(wj, 6j), i = 1,2, using (A. 11), one obtains (2.2). 

Concerning the cutoff for 61 and 62, we observe the following. (A.14) involves an 
approximation, which consists of replacing the virtual photons emitted by Z\ and Z 2 with 


real photons 71 and 72- This approximation is valid only if the masses of the virtual photons 
Ai and A2 are small compared to the mass of the produced system W. (See ref. 5, Sections 6.1 
and 6.7). By the uncertainty relations, A; « ^, * = 1,2. Hence in order for the approximation 
in (A.13) to be valid, we must have A* < W, or £ < W. Therefore, 


bi 



(A.17) 


If bi does not satisfy (A.17), the contribution to the cross section is small, and is generally 
considered negligible. See reference 5, Sections 6.1 and 6.2, and reference 12, Sections 7, 
7.1-7.3. Another consideration for the values of b; is that since we are interested in the 
effects of each nucleus acting as a single entity rather than as a collection of nucleons acting 
independently of each other, i.e., we are interested in the coherent effects of the collection of 
nucleons, we need to restrict 

bi > Ri . (A. 18) 
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So for reactions in which the Compton wavelength of the produced system is smaller than 
the nuclear radii, i.e., yfr < Ri, we can set the minimum of 6* by 

^imm == Ri • (A. 19) 

This is the case for and t + t~ pair function. But for e"*"e - pair production, the 

Compton wavelength of an electron ^ is >Ri- So we set the minimum of b{ by 

i>imm = |f- (A20) 

Hence in general we set 

^imin = max 

Note that when W = 2 m e in (A. 20), bj m , n = which is the cut off generally accepted in 
the literature for e + e“ pair production (see e.g., reference 3, Section 3). With 6j m i n defined 
by (A. 19) and (A.20), our calculations for e + e - pair production show good agreement with 
the results from the approximate formula (3.2), as is shown in Table 1. 
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Appendix B 

First we list the cross sections for the processes in (3.1): 77 — ► l + l~ , 77 — ► s 
V + V~, and 77 — + H°. (See ref. 2, eqs. (14)-(17), pp. 159, 160; and ref. 13, eqs. 
p. 95.) 


< 7(77 


i + n 


47ra 2 

W 2 


2(1 + W - lyf) ^ 


- (i + yi)\fi -ifl , 



< 7(77 


+ o- 


S ' 5 


27TQ 2 


( 1 (1 - 5 *.) ln + \/f 


< 7(77 




87ra 2 T 1 

If 2 " I £ 


1 + 7 ^ + 3 1 
4 


0 


A - 3£„(1 - 2t v ) In 



where 

_ 4 m 2 _ _ 4 m 2 

^ — jy2 ’ — jy2 ’ 


m 2 , 

tv — y^2 > A = \/l — 4£ v . 

IF is the total energy of the two photons in the center of momentum frame. 

o _2 

< 7(77 — H°) = r <5(JF 2 - mi) 

where T can be written as (ref. 12, eq. (10), p. 95) 

_ a 2 G F m z H 2 
" SirtyS 11 ’ 

and / in turn has the form (ref. 13, eqn. (11), p. 95) 


1 = S ? 2 7 <? + 5Z ^ 

q 1 


s , 77 ^ 

(10), (11), 


(B.la) 



(B.lb) 

(B.lc) 

(B.2a) 

(B.2b) 

(B.3) 

(B.4) 

(B.5a) 
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(B.5b) 


where for A > | 


and for A < | 


Iq = 3 [2A 9 + Ag(4Ag - l)/(Ag)] , 
II = 2A, + A/ (4 A/ - l)/(Af), 


I w = 3Ap^(l - 2X w )f(X w ) - 3 A^ - -, 




(B.5c) 

(B.5d) 


(B. 6 a) 


(B. 6 b) 


’’ ± = 5 ± Vi- A 

The subscripts 9 , /, and W stand for quark, lepton and W-boson, respectively. 


(B. 6 c) 


A i = ~4-, for i =■ q,l, W, (B.7) 

m H 

and m.j are the rest masses of the coresponding particles, mjj is the rest mass of H °. 

In the following, we give the derivations of the cross sections for the processes 


77 — + s + s , 


(B. 8 ) 


i + r. 


We also give a derivation of the relationship between a and T for the process 


(B.9) 


77 -> H l 


(B.10) 


Derivation for 77 — > 5 "*"3 . The lagrangian for the system, including the EM interac- 
tion, can be written as 

^EM = ~ + ieA ^ - ieA ^ <P ~ m 2 ct> + <j) (B. 11) 
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in which $ denotes a scalar field operator, denotes the photon field, p = 0, 1,2,3. We 
use the convention that repeated indices are summed over, so that for example, 


A , = A 0 A° + A\A l + A 2 A 2 + A 3 A 3 
= A 0 A 0 — A\A\ — A 2 A 2 — A 3 A 3 . 


(B.12a) 

(B.12b) 


This lagrangian can be separated into a free part, and an interaction part, so that 

Ant = ie ($ + A*-^ - A^ 

+ e 2 A»An<j> + <t>. 


(B.13) 


To simplify notations, in this derivation, we are using m instead of m 3 to denote the mass 
of the scalar particle. The 5-matrix element that contributes to (B.8) can be written in the 
form 

<p_;p+|S (1) + S (2) | k 1 ,t 1] k 2 ,t 2 > (B.14) 


in which p± denotes the momenta of s^, fcj, ej are the momenta and polarization vectors of 
the photons, i = 1,2, and S^ 2 ) is defined by 

S (2) = ^ T j Ant(*l)4nt M dS d*x 2 , (B.15) 

where T denotes the time-ordering operator. Contribution from S ^ can be represented by 
the diagrams of Figures B.la and B.lb. 

Using (B.13) and (B.15), and standard techniques of field theory, one obtains 


<p_; p + |S (2 )|fci,ei;fc 2 ,e2 > 

= ie 2 (27r) 4 (2p_ 0 • 2 p+ 0 • 2fcio • 2k 2 o)~^ O -2 

+ P— n)t 2 {p+v — kv) "b P—p) € \iP+v ~ ^v) 

k 2 -m 2 + it k' 2 -m 2 + it 

x «5 4 (p_ + p+ - ki - k 2 ), (B.16) 


where k = p_ — k\ = k 2 — p+, k' = p_ — k 2 = k\ — p+, S7 is the normalization volume, and 
t denotes an infinitessimal quantity. 
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Likewise is defined by 


and 


S (1) = i T J C M (x) d*k, 


(B.17) 


<P-;P+l^^|fci,ei;A: 2 ,C 2 > = ze 2 (27r) 4 (2p_ 0 • 2p +0 • 2fcio • 2 ^ 20 ) ly/2 
x ft -2 2ei • €2 <5 4 (P- + p+ - ki - k 2 ). (B.18) 


The diagram representing this matrix element is shown in figure B.2. 

The total cross section is obtained by squaring (B.14), averaging over photon polarizations 
e\ and € 2 , integrating over phase space, and finally dividing by the photon flux. Hence we 
have 

I 2 d 3 p-d 3 p + 


'=/i£ 

«1»«2 


<P-\P+ 


S (1) + 5 (2) 


^1, ei; * 2 , e 2 > 


(27r) f 


m2 ^ 1 

x fr x — x — , 
2c T 0 ' 


(B.19) 


in which T 0 is the normalization time. Substituting (B.16) and (B.18) into (B.19), we obtain 

42 /* ,/„\ c4 / r , \ VT 0 d 3 p - d 3 p+ ft 3 

a = e a J f{6) S (p- + p + - fcl - k + ) — 6 — 

a 2 ft 4 


= (2t)4 • 2c Q ( 2^)6 (P -° ‘ m2)V2 ^ l m Sin * X (B ' 20) 


where 


m = 1 £ 


ei ■ (k + p-)e 2 ■ {p+ - k) 


. fc 2 — m 2 + ie 
ei,f2 u 

e 2 • (A/ + p-)ei • (p+ - A/) 


+ 2ei ■ e 2 


(B.21) 

(B.22) 


k 72 — m 2 + ie 

a = (27r) 4 (2p_ 0 • 2p +0 • 2fcio • 2A:2o)"2 fi -2 , 
and 6 is the angle between p_ and the z-axis, which is chosen to be along the direction of 
k\. We shall work in the center of momentum frame of the two photons, and use the fact 
that for real photons, 
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e l • ki = €2 ■ k2 = 0. 


(B.23) 


After some algebraic manipulation, we obtain 

1 


f(9 ) = ,3 4 sin 4 0 


+ 


+ 


in which 


+ 2 - 2/? 2 sin 2 0 


e s l? - 


(1 — 0 cos 9 ) 2 (1 — /? cos 0)(1 + /? cos 9) (1 + 0 cos 0) 2 

1 1 


+ 


[1—0 cos 9 1 + 0 cos 9\ ’ 


P-o 


(B.20) can be simplified by carrying out the integration over 9, so that 


r f(9) sin 9 d9 = 4(2 - 0 2 ) + 2(0 2 + 1)(0 2 - 1)1 In 
Jo 0 


1+0 


1-0 


and hence 


a = 


e 4 a 2 f2 4 p 2 
e(lU ^0x2 7T 


(2tt) 10 • 2c 2 


4(2-/? 2 ) + 2(/? 2 + l)(/? 2 -l)-ln 


1 + 0 


0 U-/3IJ 


Using the definition of a in (B.22), we then obtain 

1 ft 4 


a = e 4 (27 t)°Q 


8n-4 


16p 4 (2tt) 10 x 2c 2 


P— o 

X 2 X 27T X 


4(2 - /? 2 )/3 


+ 2(/3 2 + l)(/3 2 — 1) In 


1 + 0 


1~0\\ 


(B.24) 


(B.25) 


(B.26) 


(B.27a) 


1 1 


(27r) 2 16p 2 0 4c 

1 + 0 


x 2tt x 


4(2 — /? 2 )/3 4- 2(/3 2 + l)(/3 2 — 1) 


x In 


1-0 


(B.27b) 


In the “natural units” in which one sets h = c = 1, this result can be written in the form 


a = 


(4tt) 2 4p 2 Q 


x 2ir 


(2 — 0)0 + (0* + l)(0 i — 1) x - In 


1 ._l 1 + P\ 


2 \l-0\ 


(B.28) 
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In terms of the variables y = W = p- 0 + p+o = 2 p- 0 , we can write 

<r = a 2 x ^ x 2 tt (1 + y)\Jl-y~ (2 - y)y In + y ^ - 1 (B.29) 

2 

in which a = |^ is the fine structure constant. This result is the same as the one obtained 
by Papageorgiu (ref. 2, eq. (15), p. 159). 

Derivation for 77 — > l + l~ . For this case the interaction lagrangian can be written in 
the form 

Ant = ~ e &( x ) 4( X M X )> (B.30) 

in which tp(x) denotes the lepton field operator. 4(x) = A^x)^ and 7 ^, p = 0, 1 , 2, 3, are 
the Dirac 7 -matrices. (See ref. 9, Appendix 2 , p. 335-361.) rp(x) = V^(x) 7 o, where ip^(x) 
is the hermitian conjugate of 4>(x). The second order term in the 5-matrix is defined by 
(B.15), with Ant( x ) defined by (B.30). The initial and final states can be denoted as 


i > = |A;i, € 1 ; fc 2 , e 2 >, 

(B.31a) 

/ > = p-,s-;p+,s+ >, 

(B.31b) 


in which we have already defined kj,ej,j = 1,2 as the photon momenta and polarization. 
p ± , s ± are the momenta and spins of l + and l~ respectively. Following the notations of 
reference 9 , appendix 2 , we can write the 5 -matrix element < /|5^ 2 ^|i > as 

< /|5^|* > = — e 2 a u(p-,s-) & v (P+’ 5 +) 

+ u(p-,s-) h -rp^~ ^YJ-~ i\ V (P+, S +) ^(P- + P+ ~ k 1- * 2 ). 
k rz — m z + ie 

where 

k = p- - ki = k2 - p+, k' = p- - k2 = k\ - p+, (B.32) 

a = (27 r) 4 ( 2 p _ 0 • 2 p +0 • 2k\Q • 2 fc 2 o) - 1 // 2 ft -2 > 
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u(p_,s_) and u(p+,s+) are the spinor wave-functions associated with l~ and l + . For 
simplicity of notations, in this derivation we use m instead of m/ to denote the mass of the 
lepton. This 5-matrix element can also be represented diagrammatically by Figures (B.la) 
and (B.lb). The total cross section is given by a formula similar to (B.19): 

E 

*l’ e 2’ 
a_ ,*+ 

Performing the sum over the photon and lepton spins, we can write 



</ 


*> 


2 <f 3 p_d 3 p+ 2° 1 

(2tt)6 2 cT 0 


(B.33) 


53 | </|‘S' (2) |i> = 53 e 4 a 2 (5 4 (p- + p+ - h - k 2 )^^{{k 2 - m 2 + ie) 2 


*-.»+ 


ei,f2 


P*Y 


x Tr|>- + m) 4x i(fi + m) fay+ - m) &(-i)(J£ + m) /ij 

+ (fc 2 — m 2 + ic) -1 ^ 72 — m 2 4- ie) -1 Tr \[$- + m) fa(k + m) fay+~ m) + m) fa] 
+ (k 2 - m 2 + ie)~ 1 (k /2 - m 2 +ie)~ 1 Tr[(jf-+ m) + m) 4i(4+ ~ rn) + wi) 4l] 
-I- ( fc 12 — m 2 + ie)~ 2 Tr [(/<_ + m) fa(k' 4- m) (p+ - m) + m) fe ] } (B.34) 

in which Tr denotes the trace operator. From (B.34), it can be seen that the sum in (B.34) 
can be naturally divided into four terms: 

T\ = Tr [(/*_ + m) i\ (ft + m) foy+ ~ m ) foi# + m ) /l] . (B.35a) 

T 2 = Tr [(jL. + m) h(ft + m) foy+ ~ m ) 4l (P + rn) fa] , (B.35b) 

T 3 = Tr [y- + m) + m) AQ4 - m) + m) 4l ] , (B.35c) 

r 4 =e Tr {{ji- +m)6(t + m) 4l(4+ ~ m) + m) h\ . (B.35d) 
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Hence we can write 


I < /| 5(2) |* > = e 4 o 2 (5 4 (p_ + p+ - ki - k 2 ) 


e 1^2 


x £{(* 2 — m 2 + ze) 2 T\ + (A : 2 — m 2 4- ze) 1 (fc /2 — m 2 + ze) 1 (T 2 + T 3 ) 


£ 1,«2 

4- (A / 2 — m 2 + ze) - 2 T 4 }. 


(B.36) 


After some straight forward though tedious mathematics, one arrives at the following: 

-^2 T l= 8 fc lo ( 1 - P cos 9) (l + + ~y~P cos 0 + p z cos z e\ - 8 m 4 , (B.37) 

4 \ ^in ^in / 


e l » c 2 


I E T 2 = i E T 3 = ^lO^ 1 - COS 2 0) [ X - 0 2 (1 - COS 2 0) 


4 ' 4 

£ 1. £ 2 e li £ 2 


(B.38) 


i 52 r 4 = 8Jfe} 0 (l + p cos 0) (1 + ^ ^ p cos 6 - P z cos 3 -8m 4 , (B.39) 

4 ei .eo ' *10 *10 ' 


z 2 2m 2 

ei ,£ 2 ' - c 10 fc 10 

where P = From (B.33) and (B.36) we obtain 


P-0 

a = e 4 (27r) 8 (2A:io) -4 £2 


_ -4/o^\8/o*...\-4n-4 J ^ 1 


(2tt ) 4 (2tt ) 6 2cT 0 

x /*/i - m 2 + ze) -2 Ti + (A ; 2 - m 2 + ze)~ 1 (A / 2 - m 2 + ze ) -1 (T 2 + T 3 ) 

J ^ « 1, £ 2 

+(Jfc ' 2 - m 2 + ie)- 2 T 4 ])<i 3 p-<i 3 p + « 4 (p- +p+ - ki - fc 2 ) 


e “ 1 f \\ Y .{^-^+^ r 2 T , 


(2n) 2 1 6fc^ 0 2 c 

+(fc 2 — m 2 + ze) - 1 (fc /2 — m 2 + ze) - 1 (T 2 + T 3 ) + (A / 2 — m 2 + ze) - 2 T 4 ]^ 


x - sin 6 d6 x 2 -tt. 
2 


(B.40) 


Now we use 


(B.41) 


k 2 -m 2 = — 2 A : 2 0 (1 — P cos 9), 
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A/ 2 — m 2 = — 2A:j 0 (l + 0 cosd), 


(B.42) 


together with (B.37), (B.38), and (B.39) to arrive at 


f ^ ^ (A : 2 - m 2 4- ie) 2 Ti sin d dd = f ^ ^ (A :' 2 - m 2 + ie) 2 T 4 sin 0 
* /0 4 «i.«2 ‘ 7 ° 4 «i,€2 

= -i/3 2 -4-8^ + 2('2 + ^Uji±4|-^ > , ( 


k* (3 \l-(3\ A:f n 1-/32’ 


(B.43) 


f ^ ^ (A : 2 — m 2 + ie) ^A ; 72 — m 2 + ie) 1 T 2 sin# dd 
J ° «1.«2 

= f ^ ^(A ; 2 — m 2 + ie) - 1 (A : 72 — m 2 4- ie) - 1 73sin0 dd 

JQ «l.«2 

= 8(1 - /J 2 ) + ^ /J 2 - 2 (2 - /J 2 )(l - /3 2 )ln|[^| . 


So from (B.43) and (B.44) we have 


f ^ ^ [(A : 2 - m 2 4- ie) 2 Ti 4- (A : 2 - m 2 4- ie) ^A : 72 - m 2 4- ie) 1 (I 2 + T 3 ) 


Cl>C2 

( z*^2 1 * \ 2 


+ (Ar* — mr + ie) ^ 4 ] sin 0 dd 


2\ii„ 1 1 + 0 1 


= -8(2 - 0*) 4- -[2 + 3(1 - ^) - (2 - /3 2 ) (1 - ^)]ln 
= - 8(1 +y) + ^-(l+y - I!-') In I -^ + 1/7 -1 , 


v " 0 V * 2 y - 1 Vy V y I 

where y = 4 jt£ = 1 — /3 2 . Putting this into (B.40), we have 


e 4 1 1 


(2tt) 2 16A:f f 


^ fc i 0 x8 -(i + y) + |(i + y-y) 


x In — — 1 x 

sfy V y J 


_e^ i_ 

(47t) 2 W 2 


x 4 tt —(1 4 - 


i+ *-t r\^ + iy- 1 


= iv2 x4,r r (1+! ' )v/1_! ' +2 (i +! '”^J ln l^ + Vv _1 IJ’ (B - 46) 

in “natural” units. This result is the same as the one obtained by Papageorgiu (ref. 2, 
eq. (14), p. 159). 
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Derivation for 77 — > aP 

s° is a neutral scalar. Using p s to denote the momentum of s°, the cross section for this 
process can be written as 

a = f\ = Y, <Ps S k h €i)k 2 ,e 2 > gflxgxi, (B.47) 

in which S denotes the ■S'-matrix. For the reverse decay process s° — * 77, the width T can 
be written in the form 

„ 1 , _ 2 d 3 fci d?k 2 o 

r = T / X, < k l,tllk 2 ,e 2 S p s > 6 Q 2 (B. 48 ) 

l ° J ei)C2 W 

From conservation of momentum, we can write 
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Therefore from (B.50) and (B.51) we now have 


<7 = 


8 c 


(2tt) 2 r x -J- 6(p so - 2fc 10 ). 
*10 


In the rest frame of s°, 


(B.53) 


2 

a = (27r) 2 r-^- 6{m 2 s - 4 k% 0 ) = — rtf(m 2 - s), (B.54) 

«10 m s 

in which s is the square of the total momentum (Aq 4- A^) 2 . For 77 —* H° in which H° is a 
neutral Higgs particle, T can be written in the form given by (B.4)-(B.7). 


28 



Appendix C. Fermion Contribution to F(H° — * 77). 

In this appendix, we derive the Fermion contribution to the decay width of the decay of a 
Higgs particle H° — ► 77. For this case the interaction lagrangian can be written as (ref. 14, 
eqs. (22.58), (22.78), pp. 676, 682) 

Ant(z) = £S( X ) + ( x )’ ( C - la ) 

where 

£ S( X ) = 4(*)tf/(*)» £ 2t( x ) = fc /tf/( x )*?( x )lM*)* (C.lb) 

tpf{x) is the fermion field operator, A^{x) the photon field operator, and tj(x) the scalar 
Higgs field operator, qj denotes the charge of the fermion, and hj the coupling between the 
Higgs scalar and the fermion. The process H° —* 77 is third order in the interaction, so that 
the relevant term in the 5-matrix is 

s<3) = I AS •^int (*l)Ant(*2)Ant(*3) <&X\ ct l X 2 ( C - 2 ) 

The initial and final states can be denoted as 

|t > = \p h >, and |/ >= |*!, ei; k 2 , e 2 >, (C.3) 


in which p ^ denotes the momentum of the Higgs scalar, kj,ej,j = 1,2, are the momenta and 
polarizations of the the photons. 

We use m j and mjj to denote the masses of the fermion and Higgs scalar. The width for 
the process is given by 




2 d?k\ <fik 2 

(2tt)6 


Q 2 


(C.4) 
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The 5-matrix element can be represented by the diagrams: 


Employing standard techniques of field theory, we find 

< * lt e i; fa,e 2 \ S (3) \p h > = -q 2 fh f (2p ho • 2 fc 10 • 2 fc 20 )- 1/2 ft“ 3/2 
x J (ftpil^ipi - fci) 2 - + ie] Tr [(pi + my) 4l(fa~ fa + rrif) 4 2 {fa ~ fa~ fa + ™f)} 

+ [(Pi “ h) 2 ~ Tr\[fa +rrif) fa + ™/) 4l(fa~ #1“ fa +™f) j 

x (p 2 — m 2 4- xe) _1 ^(pi - fci - fc 2 ) 2 - mj + iej 6 4 (fci + fc 2 - Ph )• (C.5) 

We can separate the two terms on the right hand side of (C.5) and let 5^ = 5^ + S^\ 

so that 

< * 1 , ei; fa, * 2 1 s[ 3) |p A > = -q}h f (2p ho ■ 2 k l0 • 2fc 2 o)” 1/2 fi' 3/2 * 4 (fcl + fa - p h ) 
x J (ft p\ [(pi - fci) 2 - m 2 f + *ej (P 2 - m 2 f + ie) _1 

x [(Pi ~k\~ k 2 ) 2 - m 2 + iej t{ 3 \ (C.6a) 

where 

Tj 3 ^ = T r [(jfi + THf) fi 1 + TUf) fa(fa~ fa~ fa + m /)] > (C.6b) 
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and 


< h, ei; k2, € 2 1 Ph> = ~Q 2 fhf{2ph 0 * 2Aqo • 2fc2o) _ly/,2 f2 — 3 / 2 <5 4 (A:i +k2~ Ph) 
x J d*pi [(pi - k 2 ) 2 - m 2 j + tc] [pf - + iej 

x [(Pi ~k\~ k2) 2 -m 2 + iej T^\ (C.7) 

where 

= TrJ(/<i + my) *2(/*l- #2 + *»/) ^lM" #1“ fe+ TO /)]* ( C - 8 ) 

The evaluation of the matrix elements (C.6) and (C.7) are quite similar. So we need only 
consider (C.6) in detail for illustration. By evaluating the trace in (C.6b), and using the fact 
that in the center of momentum frame of the two photons, 

Cl * ki = Cl ■ k2 = €2 • frl = e 2 • k2 = 0, (C.9) 

and also 

k\ = k\ = 0, (C.10) 

for real photons, we find 

Tp) = 4m^(4pi • ei pi • €2 — e l ‘ € 2 Pi) + 8 vrij e\ • C2 Pi ' &1 

+ 4m^(-fci • A:2 + m^)ei • C2- (C.ll) 

Now we use a standard technique of Feynman parameterization (ref. 15, Section 3.2, 
pp. 160-197). 

(pf — m 2 + ie) -1 [(pi - fci) 2 — m 2 + tej [(pi — k\ — k 2 ) 2 —m 2 + iej 
= 2 f Q dx J Q dy j^(p! - Q) 2 + pf j (C.12a) 

where 

(C.12b) 


= xk + y(tf + k%), 
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(C.12c) 


pf = -Q z -mj + kfx + (k\ + k 2 ) z y + ie 
= —Q 2 — m 2 + 2k\ ■ k 2 V + ie. 

Now (C.6) can be rewritten in the form 

< *i, ei; k 2 , ) Ph > = -<l}hf(2ph 0 ■ 2ki 0 ■ k 2 o)~ 1/2 n~ 3/2 6 i (ki + k 2 - p h ) 

x2 / ir r ' dy J cfipi [{pi - Q) 2 + pf] 3 t[ 3) . (C.13) 

FVom (C.ll) and (C.13), it is apparent that in order to evaluate (C.13) we need to compute 
the following integrals: 

*2 V = J dxd V P?Pi [(Pl ~ Q ) 2 + Pi] > 

= J dx dy d*pi pj* [(pi - Q) 2 + pf] , 

dx dy cfipi [(pi - Q) 2 + pf] . 

These integrals can be computed using the method of dimensional regularization in which 
one first computes the following integrals: 


(C.14) 

(C. 15) 
(C.16) 


/ 0 
dx dy (Ppi p^pi [(pi - Q ) 2 + pf] , 

dx dy (Ppi pj [(pi - Q) 2 + pf] , 
dx dy tPpi [(p, - Q) 2 + p?] , 


(C.17) 


(C.18) 


(C.19) 

in which n is a real number, which in the final result is allowed to approach 4. Details of 
this process is given in Appendix D. 

Prom the results in Appendix D, we find 


J d>pi 


* 

1 

So 

to 

1 

l. 

1 ^ 
i to 

[(pi - Q ) 2 + Pi 



^[ 9 '“'(p?-<3 2 ) + 4W], (C.20) 


/Al 


Pf 


.TT 2 ^ 


(Pi - Q ) 2 + pf 


3= ~ l ^r 


(C.21) 
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(pi - Q ) 2 + p?] 


7T 2 1 

3 = T 7? • 


Define 


J\ = J d*Pi [Cpi - fc i) 2 ~ m / + *«] 1 (pi “ m / + ic ) 

x [(Pi - k\ - A; 2 ) 2 - m 2 + ie j x . 


By using (C.9)-(C.12), (C.20)-(C.23), we find 

„ r l r l ~v 


J\ = -ITT 2 jf dx jf dy j^4my (p 2 - Q 2 ) 

+ 8m yQ • + 4m y k\ • fc 2 + m 2 )] • e 2 . 


Using the definitions of p\ and Q in (C.ll) and (C.12), we can simplify (C.24) to 


1 rl-x 


J\ = — i7r 2 2m t e\ • e 2 / dx 

JO JO 


4y(l - x - y) - 1 
y(l - x - y) - A -hie’ 


m3- 


where X = ^ = -f. 

The integral in (C.25) can be done after some changes of integration variables and applying 
some techniques in complex analysis. The result is 




(C.26) 


(C.27a) 


(C.27b) 
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FYom (C.5), (C.6), (C.23) and (C.26) we have therefore 

< * 1 , ei; * 2 , 62 I S^Iph > = -q) h f (2p ho ■ 2*io • 2fc 2 o)' 1/2 « _3/2 J\S 4 (h +k 2 - p h ) 

= “9/ h f (2p/io ' 2fcio ■ 2fc 2 o) -1/2 0 _3/2 5 4 (*i + fc 2 - p h ) 
x i7r 2 x 2m f x ej • £2 Jo) . (C.28) 


It is straight forward to check that 


< *1, ei; * 2 . € 2 |5^ 3 ) |j>/j > = < fci,ei;fc 2 ,e 2 |5{ 3) |p h > . (C.29) 


Therefore using (C.3)-(C.8) and (C.28), 

r “ jr \ I £ 4«/ ^ (2p to ■ 2* 10 • 2fc2o) _1 ST 3 * 4 (ki + *2 - Pft) 


e l> e 2 
.4 „ ,_2 


Now 


Therefore 


’T * im t X (‘1 ' «) 2 W>| J X ^ (2lr )6 ' 


Y. («1 • C2) 2 = 2, 2fcio = 2*20 = Pfto = > 

ei.«2 


q^fh'k 7r 4 m^ r 

r = ToTTo 2 -V 2 M>l / (*10 + *20 - P/J 

(27r) iU mir J 


(C.30) 


(0.31) 


* H 
A u2 


(C.32) 


“ leT^ «/ k / m » A i* 1 - 

is related to the Fermi coupling constant Gp by h 2 = m 2 x 2\/2Gp (ref. 14, Section 
22.2, eqs. (22.58), (22.70), and (22.83), pp. 676, 679, 684) 


r = 16^ «/ 2v/5Gf 1 J ” 1 

_ , 9/|4 a 2 Cp 3 . 2 i t i 2 

- '7 1 m ff A l J »l ■ 


(C.33) 


We note that is the charge of the fermion in units of the electron charge. Our result 
agrees with that in the literature. (See ref. 13, eqs. (10) and (11), p. 95.) We note that the 
sign of the imaginary part of J 0 in (C.27b) is opposite to that of reference 12, equation (11). 
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However, since only | J 0 \ 2 enters into quantities of physical interest, such as T and a, therefore 
this difference in sign of the imaginary part of J Q is not significant. 
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Appendix D. Evaluation of certain integrals. 


In this appendix, we outline the procedures involved in evaluating the integrals in (C.14)- 
(C.19). 

The following integrals can be evaluated by standard methods of calculus. 

roo 


f “ r _ c , x M -2t+m+l 


(D.l) 


provided m is even, m > 0, i is an integer > 0, l > y + 1, and the coefficients C m i are 
defined by 


c m ,e = 


(l 


- m - m 


1)1 f m — 1 771 — 3 3 1 


(f-1)! 


( 


l.l) 

2 2 / 


\2i—m—2 


x / (cos oy 

JO 

If £ is a half-integer, (D.l) still applies with 

/“§ 

x / (cos oy 

Jo 


(D.2) 


\2l-m-2 


(D.3) 


Using methods of complex analysis, we can show that the same formula (D.l) applies if 
M is replaced by iM in (D.l). 

J <Pp{p 2 + M 2 )~ q = J dp 0 J \p\ n ~ 2 d\p\ (p 2 - \p\ 2 + M 2 )~ a (D.4) 

where denotes the surface area of the n-dimensional unit sphere. Now using (D.l), we 
find 


/ <fp( P 2 + m 2 )-“ = r dp„(-irc^ 2A (-p 2 - M 2 r a+ ° il 

J J —00 


n-1 


= 2(— l)“3T*C n -2,o C 0a _ n-1 M~ 2Q+n . 


(D.5) 
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Prom (D.5), it is straight forward to compute 


J <f*p (p 2 + 2p-Q + M 2 ^J = j cPp (p + Q ) 2 + M 2 - Q 2 j 

= 2(-l )*? C n —2,a C 0 , a .n=l «i -2 (m 2 - 0 2 )'° +! 


(D.6) 


Now we can also evaluate 


,-<* -1 


J dTp (p 2 + 2p ■ Q + M 2 ) 


-a+1 


/<Cpp,.(p 2 + 2P « + " 2 )’“ = 2(a-l) <X?P 

= c-D ^v-r** 

xC 0a . 0 j 1 nW 2 (M 2 -<3 :! )"“ +? (? J ,. (D.7) 


Using similar techniques we can evaluate 


J <rp PuPv (p 2 + 2p • Q + M 2 ) " = (-1)*? (a - 2 - 2) (a - l)- 1 (a - 2)" 1 

X c„_ 2 , q „ 2 c 0a _^ n<,°2 2 (m 2 - « 2 )'“ +1+5 

X \ 9itu + ( M 2 - Q 2 ) (-a + 1 + |) 


x (-2 Q m Q„)] , 


(D.8) 


and 


/ 


(Tp p 2 (p 2 + 2p • Q + M 2 ) “ = (-1)^ (a - 2 - £) (a - lj-^a - 2)" 1 

X C n _2,a-2C'o ia _2|2 «?>, ( M 2 - <3 2 )'° +1+5 
n + 2(a-l-|)(M 2 -<3 2 )' 1 <3 2 |, (D.9) 


in which p^j, is the metric tensor 


9oo = “511 = ~922 = “533 = U 


(D.10) 


and all g^ v with p ^ u are 0. 
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We also note the following: 


roo u m -Id f°° u m ^ 

Jo (u 2 + M 2 ) a U ~a-ldM 2 J 0 ( u 2 + M 2)«~1 U 

Therefore, by using (D.l), we find 

C m , a = — 1 2 ^ — *) ^C m ,a—1‘ 

We can now use these results to evaluate 

J (Ppi (^PiPi - Pi9 ,i ‘') [(Pi “ Q ) 2 + Pi] 

= (-1) 2 * 3 (a - 2 - £) (a - 1 )-\a - 2)- 1 CW_2 t a-2<7 0f a-^ n i-2 


sTi 4 - n) (p? + Q 2 ) (p?) 0+5 - 4(4 - n)Q»Q v (p?) 


k -° + s 


Using (D.12) in (D.13), we find 


j <fpi (4pfpJ - p?<r) [(p, - <?) 2 + pf] ° 

= (-1) 5 ? (a - 1 - f)' 1 C n - 2 , a C 0 „.n^ fl<,°2 2 

x - n) (p? + Q 2 ) (p?)'“ +S - 4(4 - nXJ'V' (p?)"° +! 


Setting o: = 3, and taking the limit as n — ♦ 4, we find 

J rf"pi (4pfpf - p?<r ) [(pi - c) 2 + p?] ' ~ 3 

= (-1)5 X 2^2 0> 3 (*>1 “ <2 2 ) + 4< W] C 2,3 C 0,j 

= - 4 ^K(' ? -« 2 ) +4W ]' 


In similar fashion we find 


J 4”piPi [(Pi - <3) 2 + Pi] 


= (-1)°^2 (a - 1 - |) (a - 

» c aa _^ n<°! s (p ?)"“ +3 (-<?) 

= (-l) 0 * 1 2C„-2,« C 0o _^ i «‘°! 2 (p?)"° +5 (-O' 1 ) • 


(Dll) 


(D.12) 


(D.13) 


(D.14) 


(D.15) 


(D.16) 
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Therefore, 


/ Ai [(Pi - Q) 2 + A] 3 = -KQia c 0 3 nfi q“ 

. TV 2 Q* 

~ 1 2 pf 

Finally, 

/ ^pip? [(pi - <?) 2 + A ] = 2(-i) 3 / 2 <?2,j c^i n< 0) i 

. 1 


(D.17) 


(D 18) 
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Appendix E. Stopping Power 

Consider the reaction 

Z\Z 2 —* Z1Z2X, (E.l) 

in which X represents one or more particles produced in the process. Let Z\ be an incident 
particle, and Z2 represent a fixed target, whose density is p (number of nuclei per unit 
volume). Let a denote the cross section for the process (E.l), and E x the energy of the 
system X. If we disregard the effect due to recoil of Z2, then by the conservation of energy, 
the energy loss of Z\ is equal to E x ■ Consider a slab of the target Z2 of cross-sectional area 
A and thickness Ax. 


Figure E.l 


The number of Z2 nuclei in this slab is pAAx. The cross section for an incident particle 
Z\ to collide with a Z2, producing X is given by 

da 

Ac = pAAx —— dE x , (E. 2 ) 

dEx 

where we assume the energy of the produced system X to be between E x and E x + dE x . 
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Therefore the probability for this process is 


v Aa , da 

P(E x )dE x = 


dE x , 


(E.3) 


in which P{E X ) represents the probability density for the process. Therefore the total energy 
loss by the incident particle Z\ per unit length is given by 


dE 

dx 



dE x 


(E.4a) 

(E.4b) 


The — sign in (E.4) signifies the fact that energy is lost by Z\ in the process, so that the 
change in its energy dE is negative. See ref. 16, eq (6.4), page 741. 

For two-photon processes of the kind that we have considered 


E x — + U2, 


(E.5) 


in which again we use the “natural units” for which h = 1. The cross section is given by 
(2.2a). By switching the variables of integration from uq, u 2 to aq, E x , and using the fact 

dui\ du 2 = du q dE x , (E.6) 


which can be obtained from (E.5), (E.4b) can be written in the form 

— p f j — — - — E x F{ wj, E x — uj\) x <7<y«y(u;i, E x — u>\) (E.7a) 

dx J J E x — 

= p [ — [ — (wi +u} 2 ) F{u> i, w 2 )^ry( a; i> u 2) (E.7b) 

J wi J u 2 

(E.7b) is the same as (2.7). 
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Plots of W 2 








Curve 1 









Figure 4b. ®yy(W 2 ) for the reaction yy -» s + s 






Cross-section for y y V + V 
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Figure 4c. CyyOV^) for the reaction yy~»V + V" 





Cross-section for y y ~> H 
X io j (a delta function has been factored out) 

4.00 nrm i i i i i j r~i i r~i r i i i | r i n n i 1 1 | 1 n " l " ri " l " ' "i j i i i i i i 

m, op = top quark mass 


o 
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Figure 4d. ®yy(W2) for the reaction YY~> II 0 . 



P(b) for 20a Pb 20a Pb - 20a Pb 2aa Pb e 
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: These plots are only rough estimates. 



Cross-section for Ma Pb 20a Pb --> 2oa Pb Ma Pb H 



gure 6. Plots of the total cross section for the process 208p^208p^ _>2Q8 PbH® . 





y 



Figure A.L Electromagnetic fields generated by a charge q moving along the j?-axis. 
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Z2 


Figure A.2. Photons emitted by two colliding nuclei, viewed along the direction of motion of 

the nuclei in their center of momentum frame. 
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P (collision point of the photons 



Figure A.3. Cross-sectional view of the collision of two nuclei. 
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Figure B.lb 


Figure B.la— b. Second order Feynman diagrams for the process 77 — ► s + s“. 
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